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Abstract 

The quasistatic field around a circular hole in a two-dimensional hyperbolic 
medium is studied. As the loss parameter goes to zero, it is found that the electric 
field diverges along four lines each tangent to the hole. In this limit, the power 
dissipated by the field in the vicinity of these lines, per unit length of the line, 
goes to zero but extends further and further out so that the net power dissipated 
remains finite. Additionally the interaction between polarizable dipoles in a hy- 
perbolic medium is studied. It is shown that a dipole with small polar izability can 
dramatically influence the dipole moment of a distant polarizable dipole, if it is ap- 
propriately placed. We call this the searchlight effect, as the enhancement depends 
on the orientation of the line joining the polarizable dipoles and can be varied by 
changing the frequency. For some particular polarizabilities the enhancement can 
actually increase the further the polarizable dipoles are apart. 

Keywords: Hyperbolic Media, Indefinite Media, Quasistatics 

1 Introduction 



Interest in hyperbolic materials, in which the dielectric tensor is real but with its eigen- 
values taking different signs, has surged following the discovery of superlensing. The 
story of superlenses itself had its genesis in three pivotal papers. 

The first by Veselago [1] suggested that a slab of dielectric constant e — — 1 and 
magnetic permeability fi = — 1 could act as a lens. This is indicated by ray tracing, 
taking into account the negative refractive index of the slab. 

The second by Nicorovici, McPhedran, and Milton [2] (see also [3]) showed that in 
two-dimensional quasistatics one could have apparent point singularities appear in the 
field surrounding a coated cylinder with coating having dielectric constant e = —l + i8, in 
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the limit 5^-0. Specifically, with r s and r c denoting the shell and core radii, and with a 
dipole source at x = z outside the coated cylinder located in the annulus a 2 /r s > z > r s 
where a = r 2 s /r c it was proved in that paper that the complex potential V(z) outside the 
coated cylinder, where z = x + iy converged as S — » to the potential 



for all z > a 2 /z , where e c is the dielectric constant of the core. [The physical potential is 
Re{[V(,2) + V (z)]e~ %UJt / 2} where u is the frequency and t is the time.] Thus an apparent 
point singularity appears at the point z = a 2 / ' zq which lies outside the coated cylinder. 
In the limit 8 — > 0, the shell acts to magnify the core by a factor of r s /r c so it has the 
same response as a cylinder of radius a and dielectric constant e c (becoming invisible if 
e c = 1) but now the "image dipole" lies in the matrix. Within the radius a 2 / z$ > \z\ > r s 
it was numerically found that the potential develops enormous oscillations. This blowing 
up of the field within a localized region dependent on the position of the source, now 
called localized anomalous resonance, may be physically regarded as a localized surface 
plasmon and is responsible for a type of invisibility cloaking jH |3j, that has been the 
subject of considerable study |S1 d El El UHl HH H21 H31 US]. 

The third pivotal paper by Pendry [16J, which served as a catalyst for the field, made 
the bold claim that the Veselago lens would be a superlens, capable of focussing much 
finer than the wavelength of the radiation. The appearance of apparent point singularities 
caused by localized anomalous resonance was later found to justify this claim for fixed 
amplitude point sources [IZl HH UHl [201 [2U [22l [231 I2H ESI [26l [2TI [3] though a single 
polarizable dipole is cloaked rather than perfectly imaged when it is sufficiently close 
to the superlens [U [281 IH]> an d a small dielectric inclusion is at least partially cloaked 
[SI El I2H] • Despite all the work on this topic (the paper of Pendry has over 6,700 citations) 
it remains an open question as to whether large dielectric inclusions (which interact with 
the surface plasmons) are perfectly imaged when they are close to a superlens: work 
by Bruno and Lintner [6], would indicate they are not perfectly imaged (in the limit 
in which the loss in the lens goes zero) while work of [29] suggests that they may be 
perfectly imaged, though it is not clear if sufficiently small loss has been taken in this 
latter investigation. 

It was suggested by Pendry and Ramakrishna [30] that a stack of layers of equal 
thicknesses alternating between e = +1 and e — — 1, would have an effective dielec- 
tric constant of infinity perpendicular to the layers and zero parallel to the layers, thus 
channelling the field like a set of infinitely conducting wires in an insulating matrix with 
e = 0. If the two constituent materials have different thicknesses and/or the dielectric 
constants with unequal magnitudes but opposite signs then the effective dielectric ten- 
sor can be a hyperbolic material (one needs to add a small loss to the material with 
negative dielectric constant to justify this, both physically and mathematically) and the 
dispersion relation in such materials allows for real wavevectors with arbitrarily large 
wavenumbers, thus allowing for propagation of waves with arbitrarily small wavelength 
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[3T] . The breakthrough came with the independent recognition of Jacob, Alekseyev and 
Narimanov [32J and Salandrino and Engheta [33J that a multicoated cylinder or sphere 
with many thin coatings with dielectric constants of alternating signs would correspond 
to a hyperbolic material with radial symmetry and be capable of magnifying an image 
from the subwavelength scale to a scale where conventional imaging would work: the 
hyperlens was born. This was subsequently verified experimentally [3^1 135] . 

In this paper, our interest is in the two-dimensional quasistatic dielectric equation 



V • eW = 0, £ 



e x 



£, 



;i-2) 



In the hyperbolic medium we consider e x is real and e y = e x /c 2 with c = —ifi + rj where 
7] is a small positive parameter (/x and rj are real constants). In this medium, as rj goes 
to zero, (11. 2p formally approaches the wave equation, 



d 2 V 
dy 2 



dx 2 ' 



;i.3) 



and thus in this limit one should expect wavelike solutions in the hyperbolic medium 
(think of y as the time and fi as the wave velocity). In this paper we study how field 
singularities arise in the limit rj — > when there is a circular hole in a hyperbolic material, 
and we study how a pair of polarizable point dipoles interact in a hyperbolic medium 
when a uniform field is applied at infinity. 



2 The field around a circular hole in an anisotropic 
lossless medium 

Here we review the known solution for the field surrounding a circular hole in a two- 
dimensional anisotropic medium, when the dielectric constants of the medium are real 
and positive. This problem has been treated before by Yang and Chou [36] in the context 
of the equivalent problem of antiplane elasticity, but their solution is more general than 
we need and they gave the potential only as an integral. In the next section we will use 
analytic continuation to obtain the solution when the dielectric constants of the medium 
are complex. For this it is important to express the solution in cartesian coordinates 
rather than in stretched elliptical coordinates (which would be dependent on the dielectric 
constants of the medium and become unphysical when the dielectric constants of the 
medium become complex). 

Consider the transformation 

2 

z + — = 2w, (2.1) 

z 

which maps the circle \z\ — r in the z plane onto the slit — r < w < r on the real w axis, 
and maps a larger circle \z\ = r , with r > r to an ellipse E, which in the w = u + iv 
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plane intersects the axes at 

u = ±(r + r 2 /r )/2, v = ±(r - r 2 /r )/2. (2.2) 
The transformation (12. ip can be inverted to express z and 1/z in terms of w: 



iv — Vr- 

r 2 



z = 10 + V w 2 — r 2 , - = . (2.3) 



It is to be emphasized that the square root needs to be taken so that the branch cut is 
a straight line between w = r and w = — r: for computational purposes one can set 



Vw 2 — r 2 = (w + r) a/(w — r) / (w + r) , (2.4) 

where the square root on the right hand side is defined with the branch cut along the 
negative real axis. Let us suppose the material outside E is isotropic with real positive 
dielectric constant e x , and now let us make some observations. First consider the potential 
jw which at the surface z = r e ld equals 

j(z + r 2 /z)/2 = 7 (r e ie + r 2 e~ w /r )/2, (2.5) 

and has complex conjugate 

7 / r 2 r 2 r e ie \ 

where 7 is the complex conjugate of 7 and 

1/2 2 

1 / Tn T Z 



z r, 



ic I r 2 r 2 \ \/w 2 — r 2 ( r 2 r 2 



2 [r 2 + r 2 ) + 2 \r 2 r 2 J ' (2 ' T) 

Thus the potential Re( 7 u> — jt) vanishes on the boundary dE, while the potential Re( 7 w+ 
7~t) has no flux of displacement field across dE (because the conjugate potential Im( 7 w; + 
7~t) vanishes on dE). 

Let us take a parameter c which to begin with we assume is real with 1 > c > and 
let us make the additional stretching transformation 

x = u, y = v/c, (2.8) 

which transforms the ellipse which intersects the u and v axes at the points (12.21) to an 
ellipse which intersects the x and y axes at the points 

x = ±(r + r 2 /r )/2, y = ±(r - r 2 /r )/(2c). (2.9) 
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This final transformed ellipse will be the unit circle x 2 + y 2 = 1 if we choose 



r = l + C; r = y/1 - c 2 , (2.10) 

so that 

(r + r 2 /r ) = 2, (r - r 2 /r ) = 2c. (2.11) 

Inside this unit circle we put an isotropic medium with unit dielectric constant Eq = 1. 

After this stretching transformation the dielectric tensor in the exterior medium be- 
comes 

'e* 



Ey 



(2.12) 



where 

e y = e x /c 2 . (2.13) 

The potential Re(7u> — jt), expressed as a function of x and y will still vanish on the 
unit circle, while the potential Re^w + 7 ft) will still have no associated flux of displace- 
ment field across this boundary (assuming the dielectric tensor of the exterior medium 
is transformed to the anisotropic value (12.121) ) This trick of making an affine coordinate 
transformation to convert an isotropic matrix to an anistropic one has been used to find 
the solution for an isotropic sphere in an anisotropic medium (37] and to find the effec- 
tive dielectric tensor of assemblages of stretched confocal coated ellipsoid assemblages 
(Section 8.4 of [38J). More generally it can be used to obtain explicit solutions for the 
fields around three-dimensional ellipsoidal inclusions in a uniform applied field with an 
anisotropic core and anisotropic matrix, each with arbitrary orientation. 
Now consider a potential V(x, y) given by 

V(x, y) = Re(/3w + jt) for x 2 + y 2 > 1, 

= 5 x x + 5 y y for x 2 + y 2 <l, (2.14) 

where w — u + iv — x + icy and from (12.71) and (12.101) t is given by 



_ (1 + c 2 ){x + icy) 2cy/(x + icy) 2 + c 2 - 1 
1 — c z c z — 1 

At the boundary of the unit circle Re(7t) = Re(7w) so continuity of the potential V(x, y) 
requires 

Re[(f3 + 1 )w] = 5 x x + 5 y y. (2.16) 
Also continuity of the normal component of the displacement field requires 

d _ d 

n x E x — Re((3w + -ft) + n y e y^ Re(f3w + 7*) = n x 5 x + n y 5 y , (2.17) 
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where n = (n X} n y ) = (x/yx 2 + y 2 , yj \Jx 2 + y 2 ) is the unit outward normal to the 
boundary of the unit disk. Since the potential Re(7«; + jt) has no associated flux of 
displacement field across this boundary we have 

Odd d 
n^— Re(7t) + n y £ y-^ Re (7*) = n x £ x-^ Re(~7w) + n y £ yg^ Re(-jw), (2.18) 

so the flux continuity condition ( I2.17P reduces to 

.9 d 
n^— Re[(/3 - 7)10] + n y £ yg^ Re[(/3 - 7)10] = n x 5 x + n y 5 y . (2.19) 

When P and 7 are real, = and 7 = 7', corresponding to an applied field acting in 
the x-direction, then = and (12. 16ft and (12.191) are satisfied when 

f3' + i = 6 x , e x {p-i) = 5 x . (2.20) 

These have the solution 

Ox = — — r, P = — j— , (2.21) 

and from (12. 14j) and (I2.10p the potential outside the inclusion, 

V(x,y) = 0'x + -f'Re(t), (2.22) 

equates to 

j(e x + l)x j(l + c 2 )x 
V{x,y) = h 



1 -c 2 



7 c 
+ 



yj(x + icy) 2 + c 2 — 1 + yjx — icy) 2 + c 2 — 1 



(2.23) 

ir — ± 

while the field inside is 

V(s iy ) = te=?^f (2.24) 

£ x 

On the other hand when /3 and 7 are imaginary, /3 = i/3" and 7 = £7", corresponding 
to an applied field acting in the ^/-direction, then S x = and (12.161) and (12.191) are satisfied 
when 

- c(/3" + 7") = S v , -e y c(f3" - 7") = 5 y . (2.25) 

These have the solution 

S.= %**, f=?&±£, (2.26) 
and from (I2.14p and (12.1 Up the potential outside the inclusion, 

V(x,y) = -P"cy + T"Im(t), ( 2 - 2 7) 
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equates to 



V(x,y) 



i\s y + l)cy | j'(l + c 2 )cy 

1 — By 1 — C 2 



i'c 



\/(x + icy) 2 + c 2 — 1 — a/(x — icy) 2 + c 2 — 1 



while the field inside is 



i(c 2 - 1) 



1 £"71 



(2.28) 



(2.29) 



3 The field singularities around a circular hole in a 
hyperbolic medium 



The potential V(x,y) given by (I2T231 and ( IZ23J) . or by (1238]) and (127291 . solves the 
dielectric equations when c and take any real positive value. By analytic continuation 
these formulae also solve the dielectric equations when c and e x are complex, and in this 
case V(x,y) is given by substituting the complex values of c, e x and e y = e x /c 2 in these 
formula. The branch cuts in the square roots need to be taken so (for fixed non-zero 77) 
there are no singularities in the field outside the cylinder. Guided by (12. 4ft we choose 



\/ (x + icy) 2 + c 2 — 1 — (x + icy + VT 



yj {x — icy) 2 + c 2 — 1 



(x — icy + Vl — c 




x + icy — Vl — c 2 
x + icy + \J\ — c 2 

— icy — \J\ — c 2 



(3.1) 



(3.2) 



x — icy + Vl — c 2 

where the square roots on the right hand side of these expressions have their branch cuts 
along the negative real axis. The expression beneath the square root in ( 13. ip will be real 

when 

x + icy — Vl — c 2 x — icy — Vl — c 2 



x 



+ icy + \/T 



x 



icy + Vl 



(3.3) 



(where the bar denotes complex conjugation) which is satisfied when (x, y) lies on the 
line 

yRe[cVl-c 2 ] = xIm[Vl - c 2 ]. (3.4) 

Along this line the ratio in (13. 3p will be negative along the interval between the points 
where x + icy = ±Vl — c 2 , i.e. between the points 

(x,y) = ± 



RefcVT^c 2 " 



Im[VT 



Re[c] 



Red 



(3.5) 
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which will lie in the unit disk x 2 + y 2 < 1 if and only if 

(Re[c]) 2 - (Re[cVT^]y - (imf^T 3 ^ 2 ]) 2 > 0. (3.6) 

Numerically we have checked that the left hand side is always non-negative and zero only 
when c is purely imaginary. This shows that there are no branch cuts in the potential 
outside the circular hole when rj > 0. 

The parameters 7', 7", f3', (3", S x , 5 y , r and r are also generally complex. Let us take 
e x to be real and positive and c = —ifi + 77 where 77 is a small positive parameter. Then 

£y = e x /(-ifi + vf ~ -e x /n 2 + 2ie x r]/fi 3 (3.7) 
is close to being real and negative. Since 



yjx + icy) 2 + c 2 — 1 = y (x + fiy + ir/y + y/l — c 2 ) (x + /iy + ir/y — vT 



c 



2] 



a/(x — icy) 2 + c 2 — 1 = J {x — fiy — irjy + Vl — c 2 )(x — fiy — irjy — \/l — c 2 ), (3.8) 

we see that the potential V^(x, y) given by (I2.23H or (I2.28H develops singularities as 77 — > 
along the four characteristic lines 

x + fiy + a/1 + /U 2 = 0, x + fiy - a/1 + /i 2 = 0, 

x — fiy + a/1 + /i 2 = 0, x — fiy — a/1 + /1 2 = 0, (3.9) 

which are tangent to the unit disk touching it at the four points 



(x, j/) = (±1/V1 + ±/i/Vl + F 2 )- (3-10) 

Figure 1 shows a plot of the absolute value of the potential V(x, y) showing how kinks 
develop along the characteristic lines 

To better understand the behavior of the fields near these singularities let us consider 
the potential V(x,y) given by (12.231) in a region near the characteristic line x + fiy + 
+ /1 2 = but away from the disk and away from the three other characteristic lines. 
In this region V(x, y) takes the form 



V(x, y) = H(x, y) + G(x, y)yjx + fiy + irjy + r, (3.11) 



where 



Vl - c 2 « y/l + fj? + 



i fir) 



i(e x + l)x 7'(l + c 2 )x 7'ca/ (x - icy) 2 - r 2 
#1^ y) = 1 — + — 1 — o — + 2 — 1 

c — 1 



G(x,y) = . 3.12 

cr — 1 



-3-2-10 1 2 3 



Figure 1: Plot of the absolute value of the potential V(x,y) around the disk given by 
(I2.23P with an applied field directed along the x-axis and with parameters 7' = 1, e x = 3 
and c = 0.01 — i. The potential has kinks at the four characteristic lines given by ( 13. 9 j) 
which are also drawn. Near these lines the electric field is huge. 



Near the characteristic line the electric field blows up as 77 — > and 

dV uG uG 

dy y/x + ny + ir/y + r yjh + iris' 

where 

il^\J-2yJl+lJ? 
l + /i 

is the limit as r] — > of G(x, y) on the line x + f/,y + yj\ + fi 2 = 0, and where 

h = x + fiy+\/l + fi 2 , s — y-\ = (3.15) 



(3.13) 



G o = = -TVV2(1 + /, 2 )- 3 / 4 (3.14) 



(Thus /i measures the distance from the characteristic line x + \iy + yj\ + /x 2 = 0.) So 
the local time averaged power dissipated in this region per unit area is proportional to 



Im(£j 



dV 2 



dy 





\ 2 H 2 Im( 






h + irjs 





[i\Jh 2 + r] 2 s 2 



(3.16) 



in which ( 13. 7p has been used to estimate lm(e y ). Let us change variables from (x,y) to 
(h, s), so that dxdy = dhds. Observe that the right side of ( I3.15P is an even function of 
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h and that the integral 

\ - 7 ==— = Hh + ^h^7^ 2 )-Hv\s\), (3.17) 
Jo v 7r + r\ l s z 

when r) is very small, and ho is not too large (so the approximation of being near the 
characteristic line is still valid) has a dominant contribution of — 111(77 |s | ) - So when rj is 
very small we have 

2 

dhK -Ae x \G \ 2 r]\a{r]\s\). (3.18) 

Hence along the characteristic line the power absorption (integrated across the line), 
per unit length of the characteristic line, goes as 77 ln(r7|s|) which goes to zero as rj — > 0. 
However it extends a long way out along these characteristic lines so the total contribution 
does not tend to zero. To see this first note that the approximation (13. 13j) will clearly 
break down at large values of s, specifically when rjs is of the order of one, since then 
the right hand side of (13. 13[) becomes comparable to dH/dy. Thus a ball park estimate 
of the total absorption coming from this characteristic line is 

2 f ^ -Ae x \G \ 2 r ] \n(r ] \s\)ds = 8e x \G \ 2 . (3.19) 
Jo 

The interesting point is that this total absorption remains finite and non-zero as rj — > 0. 
This explains the discovery of Sihvola [39J that a hole in a hyperbolic medium may have 
loss even though the medium is essentially lossless. 



Im(e y ) 



dV 
dy 



4 The dipole approximation for the far field around 
a circular hole in a hyperbolic medium 

Consider the potential V(x, y) given by equation (I2.23j) corresponding to an applied field 
in the x-direction. When x + icy and x — icy are both large we can use the approximations 

c 2 _ y 

y/(x + icy) 2 + c 2 - 1 x + icy + — : — -, 

2{x + icy) 

\/ (x — icy) 2 + c 2 — 1 ~ x — icy H — -, (4.1) 

2 (x — icy) 



to obtain 



i{e x + l)x i{\ + c 2 )x 2icx 7'c 7'c 

V(x, y) « — + — 5 ■= + — — — — + 



e x — 1 1 — c 2 1 — c 2 2(x + icy) 2(x — icy) 



e x — 1 1 + c/ x 2 + c 2 y A 
, , e x + 1 1 - c \ / xa x 
e x - 1 1 + c J\ X x 2 + c 2 y< 



(4-2) 
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where a x is the polarizability 



a x = , ~ C = C(1 + C)(1 -^ , (4.3) 
i^e\ 2(e x + c) 1 ; 



e x -l 1 1+c 

which has been normalized to make the last bracketed expression in (14 . 2 [) as simple as 
possible. 

Similarly, when the applied field is in the y direction and x + icy and x — icy are both 
large the potential given by (I2.28P has the far field behavior 

j"(£ y + l)cy i'(l + c 2 )cy 2j"c 2 y 7 "c j"c 

V{x, y) « — + — = T + 



1 — e y 1 — c 2 1 — c 2 2i(x + icy) 2i(x — icy) 



II , 



7 cy + 



1 — e y 1 + c/ x 2 + c 2 y 



Ey + 1 1 — c \ ■y"c 2 y 



r.2 I „2„,2 



i - e y i + c y v ^ 2 + c 2 ?/ 2 



7"M^ + -)U-^)> (4-4) 



where a y is the normalized polarizability 



C _ C(l + C)(l - Eyj 



a v = -, r- = t~ (4.5) 

/ ^+l + l_c^ 2(1+05^) V ' 

We call (14. 2 p and (14. 3 j) the dipole approximation for the far field. Note that in a hyperbolic 
medium it does not suffice for x 2 + y 2 to be large to ensure that both x + icy and x — icy 
are large. One must also be sufficiently distant from the lines x = ±/j,y since along 
these lines either x + icy on - icy is close to zero when 77 is small. Thus for the dipole 
approximation for the far field to be valid one must be sufficiently far from the four 
characteristic lines (13.91) : this makes sense as the electric field diverges to infinity along 
these lines, whereas the dipole field only diverges on the two lines x = ±/xy as rj — > 0. 

The expressions (14. 3p and (14.51) for the polarizabilities could have been obtained more 
easily from the far field expressions for the potential outside an elliptical hole in an 
isotropic medium. When c is real and positive and before the stretching, the elliptical 
hole has axis lengths of 2 and 2c and in the (u, v) plane. With an applied field in the 
M-direction the far field in the isotropic medium with dielectric constant e x has potential 

u *CU(1-E X ) 

27r(n 2 + t; 2 ) [^ + (1-^)^]' 1 • J 

where 7rc is the area of the ellipse and c/(l + c) is the depolarization factor of the ellipse 
in the w-direction. By making the transformation x = u, y = v/c this potential gets 
mapped to 

X Z _|_ (yLyJ. 
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where the polarizability a x is given by (14.31) . It is similarly easy to derive ( 14. 5ft using the 
depolarization factor 1/(1 + c) of the ellipse in the v direction. When c = 1 and e y = e x 
(14. 3 p and ( 14. 5 p reduce to a x = a y = (1 — e x )/(l + e x ) which (within a proportionality 
factor) is the polarizability of a hole in an isotropic medium having dielectric constant 



5 The response of two interacting polarizable dipoles 
in a hyperbolic medium 

Here we study the mathematics of the interaction of two ideal polarizable dipoles in a 
hyperbolic medium. We leave open the question as to whether these ideal polarizable 
dipoles have any physical significance. Nevertheless the searchlight effect discussed here 
should motivate future studies to see, say, whether a distant very small circular disk, if 
appropriately positioned, can substantially influence the response of a large circular disk. 

By definition a polarizable dipole with rectangular symmetry located at the origin 
responds to a local field acting on it in the x direction so that the potential V(x, y) close 
to the origin has the expansion 

a x /2 a x /2 



V(x, y) ~ x + a\ — 



X + <2i — 



(5.1) 



x + icy x — icy ~ x 2 + c 2 y 2 ' 

and responds to an local field acting on it in the y direction so that the potential V(x, y) 
close to the origin has the expansion 

a y /(2ic) a y /(2ic) _ ya y 



V(x,y) &y + a 2 + 



y + a 2 - 



x 2 + c 2 y 2 



(5.2) 



x + icy x — icy 

Here we call a x and a y the polarizabilities of the polarizable dipole, a± and a 2 are con- 
stants, and c = ifi — 7] with rj being small. By taking linear combinations, the response 
to an arbitrarily oriented local field is such that the potential close to the origin has the 
expansion 



V(x, y) w >y x x + i y y + a 



x 2 + c 2 y 2 



(5.3) 



With c = — i/j, + f] and rj > small, the potential on the right hand side of ( 15. ip now 
has a local time averaged power dissipation near the characteristic line x + fiy = of 

2 



lm(e l 



dV 



dy 



» 2 \ 


a x 


2 Im( 




4\g + ir]y 


4 



v £ x\<y x 



2fi(g 2 + rfy 2 ) 2 ' ^ 

where g = x + jiy. So when rj is very small the dissipation integrated with respect to g, 
in the range — go > g > go, is approximately 



Im{e y 



dV 
dy 



dg 



2fiy 3 r/ 2 



du 



(5.5) 
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where v = g/(t]y)- Thus this power dissipation integrated across the characteristic line, 
per unit length of the characteristic line, blows up as r] — > 0. 

Now consider a uniform applied field in the x-direction acting on two polarizable 
dipoles, each with rectangular symmetry, one located at the origin and the other at the 
point (x , yo). The total field V(x, y) is the sum of the uniform field plus the two dipolar 
fields: 

Xplx + Vply , (% - X o )0 2x + (y~ yo)p2y 



V(x, y) = x + 



+ 



x 2 + c 2 y 2 (x — Xq) 2 + c 2 (y — y ) 2 
Expanding this around the origin x = y = gives 

xfkz + yPiy (x (3 2x + y o 2y ) 



(5.6) 



V(x,y) 



x 2 + c 2 y 2 Xq + c 2 yl 

x(3 2x + y^2y (2xx + 2yy )(x (3 2x + y /3 2y ) 



+x + 



x 2 + c 2 y 2 



{x 2 + c 2 y 2 ) 2 



(5.7) 



which allows us to identify the local field acting on the dipole at the origin. If (ai x , cti y ) 
are the polarizability coefficients of the dipole at the origin then from ( 15. 3 p we have 







lx 



-1 - 



2:v 



xl + c 2 yl 



+ 







2xq(xqI3 2x + y (3 2y ) 
{xl + c 2 y 2 ) 2 
2c 2 y (x (3 2x + y o 02y) 



x 2 + c 2 y 2 



{x 2 + c 2 yl) 2 



a 



(5i 



In a similar fashion, by rewriting ( 15. 6 p as 



V(x,y) = x + (x- x ) + 



[x + (x- x )]Pvb + [yo + (y- Vo)]0i y . {x - x )(3 2x + (y - y )/3 2y 



[x + {x- x )] 2 + c 2 [y + (y- y )]' 



+ 



(x - xq) 2 + c 2 (y - y )< 



(5 



and expanding this around the point (xo,yo) we obtain 
V(x,y) ; 



x - x )(3 2x + (y- y o )02y . . (x (3 lx + y o 0i y ) 

+ Xq + 



x - x ) 2 + c 2 (y - y ) 2 

(x - x )/3 lx + (y- y o )0i y 



x 2 + c 2 yl 



+{x - x ) + 



x 2 + c 2 yl 



[2(x - x )x + 2(y - y )yo]{x p lx + y Pi y ) 
(xl + c 2 yl) 2 



(5.10) 



So if (a 2x ,a 2y ) are the polarizability coefficients of the dipole at the point (x ,y ) then 
we have 



02 X 

fa, 



-1 







lx 



x 2 + c 2 yl 



+ 



0iy 



Xq + C 2 yl 



+ 



2x o (x o lx + y o ly ) 
(xl + c 2 yl) 2 
2c 2 y o (x o lx + y o ly ) 



Ot2x 



(x 2 + c 2 yl) 2 



a 2y . 



(5.11) 
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Introducing 

e = %l + c 2 y 2 Q ^x 2 - n 2 y 2 Q - %mvl, ( 5 - 12 ) 

in terms of which c 2 = (e — xfy/y 2 , the equations (I5.8P and ( 15.1 ip take the equivalent 
form 



Ax 
OCi x 

Ay 
ai y 

Ax 
a>2x 

Ay 

a 2y 



-1 



Ax 2x (xq/3 2x + y l3 2y ) 
e 



^ 2(e - xg)(x /3 2x . + ?/oAy) 
e ?/ e 2 

0ix 2x (xq/3i x + ?/oAy) 
e 

Ay 2(e - xg)(x Ax + yoAy) 
=2 



-1 



These four equations have the solution 



(5.13) 



(3 



where A is the matrix 



A 



i 



"Ax" 




"-1" 


Ay 


= A- 1 





Ax 




-1 


-Ay. 




. . 





i _ 


2x 2 


e 




i 


2x 


o(e-x 



2x y 



1 

"2.1 



2x (e-x 2 ) 



1 I 2 xq 
"7 "I 7T- 



which has determinant 



2xo|/o 



1 

a 2y 



det(A) 



[4x 2 (e - x 2 )(aix - a ly )(a 2x - a 2y ) + (e 2 - a lx a 2x ) (e 2 - a ly a 2y )} 



e*ai x ai y a 2x a 2y 
that vanishes when x^ solves the quadratic 

4z 2 (e - x 2 )(a lx - a ly )(a 2x - a 2y ) + (e 2 - a lx a 2x ) (e 2 - aiyt^y) = 0. 



(5.14) 



(5.15) 



(5.16) 



(5.17) 



We are interested in what happens to this solution when x an d yo are such that e is 
very small. Using Maple to compute the matrix inverse and taking the limit e — > we 
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(5.1* 



find that 

2xo 2 ai x ai y \2 xo 2 {a2y-a2x)+ct2ya2x\ 
a ly a2ya lx oi2 x -4:Xo 4 (a lx -a 1 y)(a2 x -a2y) 

2xo 3 Ql x ai y [2 XO 2 (ct2y—Ol2x)+ai2ya2x\ 

yo [ai y a2 y ai x a 2x -i x 4 (a lx -a ly ) (a 2x —OL2 v )] 
2x( 2 a.2 x OL2y [2a;o 2 (oii,-aii:)+ai :c oij / ] 

CL\yOl2yOI.\ x a2 x — 4X0 4 (oia;— Ct\y)(u2 x —OL2y) 

2xo 3 a2 y a2 x [2 xo 2 (ai y -a lx )+a\ x aiy\ 
yo [a\ya.2yCt\ x a.2 x — 4 xo 4 (ai x — ai y ) (pc2 x —Q.2 y )] 

Thus when a.\ x ^ a. ly and a 2x ^ ct 2y we see that in the limit e — > all components of /3 
blow up to infinity when xq is such that 



Xn 



a lx ai y a 2x o:2y 



4:(a lx - a ly )(a 2x - a 2y ) 



(5.19) 



which is in agreement with (I5.17P when one sets e = 0. Of course this equation generally 
will not have a solution for real Xq if any of the polarizabilities are complex. 

Additionally let us suppose that the polarizabilities a 2x and a 2y are very small, and 
«2x 7^ Oi 2y . Specifically let suppose that a 2y = ka 2x , where k is a fixed constant not equal 
to 1. Then if the limit a 2x — s- is taken after the limit e — > 0, ( I5.18P implies 



(3^ 



ot\ x —a.\y 
ai x ai y XQ 
(a lx -a ly )y 






(5.20) 



By contrast if a 2y = ka 2x , where A; is a fixed constant, and we take the limit a 2x — > 
directly in f!5.14p . keeping e fixed and non-zero, then we obtain 



-Mix 







(5.21) 



Thus even a polarizable dipole at (xo,yo) with very small polarizability can have a very 
large effect on the net dipole moments of the system if it is positioned close to one of 
the characteristic lines where e is small. Even though its polarizability is small the field 
exerted by this dipole on the dipole at the origin is still significant. We call this the 
searchlight effect since fi and thus the angle of the characteristic lines will depend on 
frequency, so by varying the frequency and observing the net dipole moments of the 
system one may hope to detect something about the relative location of the polarizable 
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Figure 2: Plot of the absolute value of the dipole amplitude /3i x near the line where e ~ 
with polarizabilities a\ x = 2, and a.\ y = 1, a^ x = 0.2, and a<i y = 0.1 and parameter 
c = 0.01 — i. We use rotated coordinates £ = (x + y )/y/2 and r = (x — yo)/y/2. Note 
the long-range interaction. 




Figure 3: Same as figure 2 but with polarizabilities a\ x = 2, and a\ y = 1, ot2 X = 0.1, and 
= 0.2 and parameter c = 0.01 — i. 
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dipoles, even though one of the polarizable dipoles, by itself, is difficult to detect. The 
effect is illustrated in figures 2 and 3. 

The case k = 1 when a 2y = a 2x is rather special, but very interesting. In this case, 
with e — > (I5.18P simplifies to 



/3-> 



2x 2 

2XQ 3 

yo 

2 x 2 x 2 ai x +2x 2 ai y ) 
aixOtiy 

2 xq 3 (ai x a-Ly-2 xq 2 ai_ x +2 xq 2 aiy} 
yoa ly a lx 



(5.22) 



and this result does not depend on the magnitude of a 2x . Remarkably, note that the 
magnitudes of the components of f3 increase as xq increases: recalling that when e = 0, 
yo = ±xo/(ic) we see that fii x and (3i y increase as Xq, while f3 2x and (3 2y increase as Xq 
when a\ x ^ a± y and as x^ when a\ x = a\ y . Thus the interaction increases the further 
the polarizable dipoles are apart! 



abs(/3 lx ) 




Figure 4: Plot of the absolute value of the dipole amplitude (3\ x near the line where e ~ 
with polarizabilities a\ x = 2, and a\ y = 1, a 2x = a 2y = 0.1 and parameter c = 0.01 — i. 
We use rotated coordinates £ = (xo + yo) / \/2 and r = (x — yo) / \/2. Note that the 
interaction along the line r = first becomes stronger as £ increases, then weakens. 

To shed more light on this one can, using Maple, directly compute the right hand 
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side of (I5.14p when a 2y = ct 2x to obtain 




(e 2 -a 2x a ly )y 
<X2xf 



(5.23) 




y (,e 2 -a2 x a ly )(e 2 -a lx a 2x ) 



where 



9 



e 3 ai x -2x 2 ai x e 2 - ea lx a 2x a ly + 2x 2 a lx a 2x a ly - e 4 

2 12 2 4 

+4a lx a 2x ex -Aai x a 2x x + a 2x a ly e - 4 x ct\ y a 2x e + 4 x a ly a 2 
a lx e 2 - a ly a lx a 2x - a lx a 2x e + a ly a 2x e + 2 a lx a 2x x 2 - 2 x 2 a ly a 2x . 



(5.24) 



If the limit e — > is taken in this expression we recover ( 15.22ft . On the other hand 
it is evident from ( 15. 23ft that there are resonances when e 2 equals ai x a 2x or ot\ y a 2x) and 
that the resulting expression for (3 depends crucially on the ratio of the magnitude of 
e 2 to these two quantities. Also if r\ is non-zero the interaction decreases for sufficiently 
large separations of the polarizable dipoles. On the characteristic lines x = ±fiyo we 
have e ~ —2i\ir\y\. So for large xq = ±fiyo ( 15.23ft implies that, for example, 



which goes to zero as l/y^ as y — > oo. Figure 4 shows how, as the separation increases, 
the interaction along the characteristic line first increases, then decreases. 
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